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rate for this process being calculated by a method of moments. Our calculated spin decoherence 
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serious problem in developing spin-based semiconductor quantum computer architectures. 
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I. INTRODUCTION 

Electron and nuclear spins in semiconductors are 
promising qubit candidates for quantum computation be- 
cause their intrinsic quantum two level nature together 
with existing semiconductor microelectronics technology 
can potentially satisfy the strict control and scalability 
requirements of a quantum computer (QC). Hence elec- 
tron spins in quantum dots (QD) 1 and donor impurities 2 
as well as nuclear spins in semiconductors 3 have been 
suggested as candidate building blocks for feasible QC 
architectures— However, to build such a device major 
advances in single spin manipulations are needed, and 
for this purpose realistic calculations of semiconductor 
spin dynamics are essential to guide the experimental 
effort currently taking place. A question of particular 
importance is whether a localized spin will remain unaf- 
fected by the many interactions invariably present in a 
semiconductor environment during a time interval long 
enough for fault tolerant quantum computation (equiv- 
alent to 10 4 — 10 6 quantum gating times. 5 ) In a recent 
paper— we showed that spin coherence of bound electronic 
states in semiconductors is limited by spin-spin interac- 
tions at low temperatures. When this interaction is be- 
tween the qubits themselves, it can in principle be incor- 
porated into the QC Hamiltonian, although this will lead 
to more complicated gate sequences. In particular we are 
not aware of any theoretical QC work specifically work- 
ing out such gate sequences incorporating inter-qubit in- 
teractions. Therefore it is instructive to analyze the er- 
ror introduced by ignoring some of these interactions, 
as we did in the case of dipolar coupled spin qubits. 6 
The presence of many non-resonant spins in the system, 
such as lattice nuclei, also leads to phase fluctuation of 
the spin qubit, an effect which is hard to control. This 
has been denoted Spectral Diffusion (SD) since the qubit 
Zeeman frequency will diffuse through the spin resonance 
line. Spectral diffusion specifically refers to fluctuations 
in the Zeeman frequency ^B g = g^iBB c g/h (where g is 
the effective g- factor, fis — eh/2mc the Bohr magneton, 



and B c g the effective local magnetic field) of the electron 
due to external effects arising from the semiconductor 
environment. Note that such fluctuations could arise ei- 
ther because the effective magnetic field B e s is changing 
dynamically or because the electron ^-factor is varying. 
There are many physical processes leading to spectral 
diffusion, and here we are specifically interested in the 
limiting processes at the lowest temperatures. The phys- 
ical process of interest to us is therefore dipolar nuclear 
fluctuations. 

The first order of magnitude estimate of this effect was 
applied to Si:P donor electrons^ while our recent paper 6 
used the same methods to estimate the SD rate in GaAs 
QDs. However, these estimates assumed a priori that nu- 
clear pairs flip-flopped slowly (and hence the echo decay 
was r~j exp — t 3 ) with a rate given by an unjustified phe- 
nomenological equation [Eq. (15) in Ref . flTI and Eq. (8) 
in Ref. ||| . Here we propose a new description for this de- 
coherence mechanism, arising from dipolar fluctuations 
of the lattice nuclear spins, affecting the qubit Zeeman 
frequency through hyperfine coupling. Even though we 
still treat each nuclear pair as an independent Markovian 
random variable (an approximation which seems reason- 
able for temperatures well above nuclear dipolar order- 
ing, happening on the nanokelvin scale), our new theory 
describes fast and slow flip-flops on the same footing in- 
corporating motional narrowing effects previously absent 
in former treatments (which happens when the fluctua- 
tion is so fast that SD is suppressed) . We also derive mi- 
croscopic expressions for these flip-flop rates, leading to a 
more refined calculation of nuclear SD for GaAs QDs and 
Si:P donor electrons, together with the first treatment of 
this effect for a 31 P donor nucleus. 

In the case of localized spins precessing in a magnetic 
field B, knowledge of three phenomenological parameters 
is sufficient to describe the spin 1/2 dynamics: The gy- 
romagnetic ratio 7 which determines the precessing fre- 
quency (or equivalently the g factor, g — 2mcj/e with e 
the electronic charge, m the bare electron mass and c the 
speed of light) , the longitudinal relaxation time or spin- 
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flip time T\ , and finally the transverse relaxation time or 
dephasing time T 2 , which is often denoted coherence time 
since it sets the time scale of coherent superpositions be- 
tween states along the B field direction^ However, elec- 
tron spins in a solid state environment often have differ- 
ent precession frequencies, either due to hypcrfinc fields 
from nearby nuclear spins or from unequal gyromagnetic 
ratios (arising, for example, from varying carrier effective 
mass). Therefore the transverse magnetization of a spin 
ensemble will decay in a time scale T 2 which is in general 
much shorter than the single spin dephasing time T 2 . The 
latter time scale can be measured using a tt/2 — tt spin 
echo sequenced The time it takes for this echo to decay to 
1 /e of its initial value conveniently defines our single spin 
coherence time and has been historically called Tm (spin 
memory time)^ since the echo envelope usually does not 
decay exponentially as predicted by the Bloch equations 
from which T 2 was first defined^ (In appendix B we show 
that measuring a tt/2 — tt echo is equivalent to measuring 
the modulus squared of a single spin off diagonal density 
matrix element). 

It has been known for a very long time that SD caused 
by nearby non-resonant spins is usually the dominant 
echo decay mechanism— However, all former SD theo- 
ries assumed a single relaxation rate for the non-resonant 
spins, an approximation perfectly suitable for "Ti sam- 
ples" , whereby these spins change their states through 
spin-flips only. The theories of B. Herzog and E.L. 
HahnjiS and later J.R. Klauder and P.W. Anderson, 11 
described the central spin Zeeman frequency as a ran- 
dom variable evolving in time according to Gaussian and 
Lorentzian conditional probabilities respectively. These 
assumptions lead to a tt/2 — tt echo decay of the form 
exp (-2Tf 1 J 2 r 3 /3) and exp (-Tf 1 St 2 ) respectively as 
long as t <C T\, with 2r being the time interval be- 
tween the first pulse and the echo. Both 8 (t — > oo 
linewidth for the conditional probabilities) and T\ (spin- 
flip time of non-resonant spins) are parameters that can 
in principle be calculated from the interactions. If the 
condition r ^> T\ is satisfied, one obtains exp (— T\S 2 t) 
and exp {—St) respectively, characterizing the motional 
narrowing regime. Interestingly, Gaussian SD correctly 
describes motional narrowing since T\ appears in the 
numerator, but Lorentzian SD does not, the decay be- 
ing independent of T\. Motivated by this inadequacy 
of the Lorentzian theory, G.M. Zhidomirov and K.M. 
Salikhovsi^ proposed a many parameter model, which 
treated the number of flips of a spin i during a time inter- 
val t as a Poisson random variable parametrized by t/Ti, 
the frequency change on the central spin being A j . Their 
theory obtained the correct motional narrowing limit and 
agreed with experiment in dilute T\ samples where the 
non-resonant spins are randomly distributed. Our prob- 
lem, however, is in a completely different regime. Here 
SD is caused by Si and GaAs lattice nuclei, which have 
Xi of the order of hours and hence the relevant time scale 
is given by the dipolar interaction, which varies substan- 
tially depending on the specific pair flip-flopping (such 



a system is denominated a "T 2 sample", since T 2 <C Ti 
for the spins that create the SD effect). We generalize 
the latter theorjiiS to many relaxation rates X," 1 , each 
corresponding to a pair n,m of nuclear spins. We also 
present a microscopic theory to calculate these flip-flop 
rates. Altogether this approach is, to our knowledge, the 
first systematic attempt to describe SD in T 2 samples. 
We also show that our theory reduces to the earlier sim- 
ple approximations of Refsltallllin the appropriate limits. 

The phosphorus donor impurity in silicon is the text- 
book example of a localized electron spin in a semicon- 
ductor. It has been extensively studied experimentally 
using Electron Spin Resonance (ESR)i2iiiii£ and suc- 
cessful theories for its gyromagnetic ratio and T\ were 
developed. 15 ' 16 However Tm for P in Si remained unex- 
plained, even though it was measured thirty years ago— 
Our model leads to a Xm two times longer than the mea- 
sured value (this agreement should be considered reason- 
able since existing theories for T\ are also off by a factor of 
two— ai 6 -) Our theory predicts a smooth transition of the 
echo envelope from Gaussian SD to motional narrowing 
behavior, this transition being well described by a cor- 
relation function. If the echo decay of 31 P donor nuclei 
is measured, we predict an echo envelope purely expo- 
nential, well into this motional narrowing regime, with 
Tm = 0.60 s. An important point discussed here is how 
SD is rapidly suppressed by reducing the amount of nu- 
clear magnetic moments in the lattice. Our calculations 
show that isotopic purification of Si (exchanging spin-1/2 
29 Si nuclei by spin-0 28 Si) may lead to coherence times as 
long as 100 ms for P impurities in Si, a result supported 
by recent experiments— Unfortunately, Ga and As nu- 
clei have no stable spin-0 isotopes, hence it seems that 
the only way to increase spin coherence in these materials 
is to suppress flip-flop events by nuclear polarization, as 
can, for example, be done by applying a strong external 
magnetic field or by using the Overhauser effect. Fur- 
thermore, we recently reported the first Tm calculations 
for a GaAs quantum dot£ This was particularly impor- 
tant since Tm has never been measured in this system, 
and a realistic assessment of the feasibility of a quantum 
dot quantum computer was needed. The detailed cal- 
culation presented here confirms our previous estimates. 
Hence the present paper together with other recent cal- 
culations of g factor— and Xi— available in the quantum 
dot literature provides a general picture for electron spin 
dynamics in these heterostructures. 

It is instructive to clarify the relationship between our 
results and recently published theories^ on related is- 
sues. In Ref. the authors considered a Hamiltonian 
which contained only hyperfine couplings between a sin- 
gle electron and the lattice nuclei, discarding the es- 
sential ingredient of the spectral diffusion effect, which 
is the dipolar interaction between nuclei. Hence their 
mechanism is based on flip-flops between electron and 
nuclear spins. But when a B field is applied electron- 
nuclear flip-flops are forbidden by energy conservation, 
since the nuclear Zeeman energy is 10 3 times smaller than 
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the electronic Zeeman splitting. Therefore their mecha- 
nism is only relevant at low B fields, when the hyper- 
fine coupling is of the same magnitude or greater than 
the electronic Zeeman energy, leading to the condition 
B < ^7/|*(0)| 2 < 100 G where 7/ is the nuclear gy- 
romagnetic ratio and (^(O)! 2 = 10 22 — 10 25 cto~ 3 is the 
electron's probability density on a nucleus. The theory 
presented here is valid in the opposite limit, B 100 G. 

This paper is organized in two parts: General theory 
and applications. In the first part (sec. II) we describe 
our theory of spectral diffusion due to a dipolar coupled 
spin system. This general theory can be easily applied 
to other spin resonance experiments, such as three pulse 
echoes. In the next part (sec. Ill) we give numerical 
results for three particular cases, and discuss their im- 
plications for the current experimental effort in semicon- 
ductor spin quantum computation. We conclude in sec. 
IV with a summary and some general comments. 



II. GENERAL THEORY 
A. Stochastic theory for the nuclear bath 

Our problem is to describe the dynamics of a carrier 
spin S coupled to a lattice of nuclear spins I„. The 
total Hamiltonian can be separated into three parts: 
H = Hs + 7~Lsi + T~ti, 



H s = j s BS z , 



Ti-si = ^ A n I nz S z 



(1) 



(2) 



til — —JlB/]Inz ~ 4 bnmlnzlmz 
n n < m 



(3) 



A n — 



,„- ; /Hy|*(R„)| 2 

d 3^( r) | 2 ^_ R «l 2 - 3 [( r - R ")-] 2 



r-RJ 



1 y ;*(Rn): 2 



1 - 3cos 2 n l 



1 2 l-3cos 2 6>„ m 
Km = -jiih jr 3 • 



4) 



(5) 



(6) 



Here the Hamiltonians are divided by h to simplify the 
notation; 75 and 77 are gyromagnetic ratios, A n the 



coupling with a nucleus located at position R„, R nm 
the relative vector between two nuclei and 6 nm the an- 
gle between this vector and the B field direction. The 
electron-nucleus coupling A n includes a Hyperfine term 
and a residual dipolar interaction. The Hyperfine term 
comes from the singularity of the integral [Eq. , which 
is removed by integrating over the angular coordinates 
first. Here we will assume this dipolar term is only ef- 
fective for |R„| > ro, which is a proper length scale for 
the electron's wavefunction (9 is the step function, while 
n the angle between R„ and the B field). The nuclear 
spins are in constant turmoil due to their mutual dipolar 
interaction. To see how this affects the spin S we ap- 
proximate each nuclear spin operator 2I nz by a classical 
random variable cr n (t) = ±1, which is valid only if the 
nuclei have spin 1/2 (this description is still accurate for 
spin / > 1/2, as long as the nuclei are in the slow SD 
regime, see below). Hence the Zeeman frequency, u) z , of 
the electron spin becomes 



On the other hand, the evolution of the nuclei is strongly 
affected by the field produced by the central spin S. This 
effect is treated by assuming the nuclei evolve according 
to the effective Hamiltonian 



H'j = Hi + - A n I n 



(8) 



which conserves total spin in the z direction. Therefore 
when any I nz flips, a corresponding I mz must flop in the 
opposite direction. These flip-flop events show that we 
can not treat the random variables a' n as independent 
of each other. Rather, we will treat pairs of spins as 
independent random variables. Hence Eq. Q becomes 



with A, 



= X! A„ m cr„ m (i) + const, 



\A n - A m \/2, and a n 



(9) 



: ±1 random vari- 
ables uncorrelated with each other. We further make 
the Markovian assumption that the probability that a nm 
changes sign during a time interval t is given by t/T nm , 
independent of past values of a nm (this Markovian ap- 
proximation is reasonable in the absence of any contrary 
evidence about the stochastic fluctuations of the nuclear 
spins). Hence the number of flip-flops N(t) is a Pois- 
son random variable with parameter t/T nmi and we may 
write 



with N(t) having the distribution 



P(N(t) = k) 



t 



1 



w I 7 

^ • \ J n n 



exp 



(10) 



(11) 



In the next section we show how to calculate the flip-flop 
rate T~l. 
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We now proceed to the derivation of the spin echo de- 
cay. The complex in-plane magnetization 



broadened line, not an echo peak>i2 The slow SD limit 
1 x implies 



nrn ^ 



v(t) = (S x )+i(S y ), 



(12) 



can be calculated for any spin echo sequence by taking 
the averag o 10 ! 11 

v(t) = ^cxp (i J s(t')Lu z (t')dt'^j ^ . (13) 

For the ir/2 — it echo considered here the echo function 
s(t) = 1 for < t < t and s(t) = — 1 for r < t. Therefore 



){t) = Y[ V nm{t), 



(14) 



exp 



+ 



cos [A nm (t- 2r)] 



1 



A T 
£>(A~ 2 T~ 2 ) 



sin (A nm t) 



(20) 



After performing the product over these slow nuclear 
pairs, Eq. (|14|) will have a peak at t w 2r. This is 
because if t ^ 2r, the argument of the product will be 
zero for some pair n, m making the whole product vanish. 
Assuming the echo peak occurs exactly at 2r, we get 



(2r) « exp <^ - 



1 



Tn 



2t- 



sin (2A„ m r) 



(21) 



cos 



A 



nm / ^ 





(t'X-i^'W 



(15) 



where we take an average over <r nm (0) = ±1 with prob- 
ability 1/2 (this average applies to an ensemble of spins; 
however, we show in Appendix B that the effect of the 
echo is precisely to remove this average, making v(t) ex- 
actly equal to a single spin in-plane magnetization). In 
Appendix A we calculate this average, and the result is^ 



Vnm (t) 
<*>(*) 



6(T-t)v£l(t)+6(t-T)vW(t), (16) 



exp 



t 



T 

± n 

COSh (R nm t) 



1 



■ sinh (Rnmt) 



(17) 



nm \ L J 



Rnm ex P 



t 



T 2 



■ COSh (R nm t) 



Rr, 



■ sinh (Rnmt) 



-A 2 m cosh [R nm (t-2r)} 



(18) 



with R 2 im = T~ 7 2 — A 2 m . We distinguish two limits in the 
above expressions: For nuclear pairs causing fast spectral 
diffusion, consider 3> A nm . In that case we have 



,(E) 



(*) 



(t) 



exp 



^A 2 T t 

2 nm- L nm b 



-l-C(A nm X 1 nm ), 



(19) 



which shows that fast flip-flopping nuclei do not even 
form an echo. Rather, they contribute to a "motional 
narrowed" signal, 8 which decays exponentially in a time 
scale that lengthens as increases. The idea is that 
the echo sequence can not refocus the spins if their Zee- 
man frequency changes appreciably over r. For large 
the nuclear pair contributes to a homogeneously 



If in addition to this slow spectral diffusion regime, we 
look at the limit A nm T <§; 1, the result is 



(2r) « exp 



1A? 



6 T n 



-(2rf 



(22) 



which is very similar to an expression derived previously 
by us& using Gaussian spectral diffusion as a starting 
point. Apart from a factor of 2 these expressions differ 
by a term which takes into account the small broadening 
due to spin-flip processes. Since this term was introduced 
heuristically and only changes numerical values by a neg- 
ligible amount we will not include it here. 

Notice that nothing should be concluded about the 
qualitative decay of v(2t) before performing the prod- 
uct over pairs in Eq. (|14fl . For example, in the particu- 
lar case of SD caused by dilute paramagnetic impurities 



(A, 



with r being the impurity-electron distance) 



it can be shown that v(2t) ~ exp (-at — br 2 )— after 
calculating Eq. I|14l) and taking a spatial average. In- 
deed exactly this behavior was seen in the electron spin 
echo decay of Si:P when the P concentration was high 
enough such that spectral diffusion due to nearby non- 
resonant electrons was dominant. 17 Here we deal with 
even more complicated expressions for A nm [our Eqs. 
(|57|) . (|rJ2|) . and JSJ}]. Moreover it is often the case that 
T,^ 1 ~ A nm and use of the limits (|T§|) , l(2T|) becomes un- 
justified. Hence in our calculations below we perform the 
product 114fl using the exact expressions 1|18|) . 

Finally, it is easy to calculate the correlation function 
[see Eq. ijAHl ] 



Xc(t) = lim 



M*Mt')> - M*)>Mf)> 



f-rt (L0 2 (t')) - (L0(t')) 2 

2t 



Pn 



Pnm ^ X P ( rp 

, \ J- nm 

Vr'A 



T -1 A 2 

nm nm 



-1 A2 

y ' 



(23) 
(24) 

(25) 
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The averages in Eq. lf^|) assume a nm (t = 0) = +1 for 
all pairs n,m. Then Xc(^c) = e_1 defines the correla- 
tion time for the effect of the nuclear bath on the spin 
S. While the two parameter Gaussian SD theory has a 
simple exponential correlator \c — exp (—t/Ti), Eq. I|24|) 
in general does not decay exponentially. This correlator 
makes clear the difference between our many parameter 
(Tnm, A nm ) stochastic theory and these simpler theories. 
If t t c , we will have Xc(t) ~ 1, each a nm (t) performing 
less than one flip-flop on average. Below we show nu- 
merically that this regime leads to ve(2t) ~ exp (— r 3 ). 
Therefore, when our calculated Tm <C t c , we predict 
a Gaussian SD decay for the echo. However, when 
Tm 3> t c , Xc(Tm) <C 1 indicating that most a nm have 
already performed many flip-flops. In this case we ob- 
tain ve(2t) ~ exp (— r), characterizing a motional nar- 
rowing regime since the nuclear bath dynamics may be 
considered fast. Therefore the importance of our corre- 
lation function is to describe the transition to motional 
narrowing induced by all nuclear pairs, not just a single 
one satisfying T m \ > A nm . 

It is important to discuss the case of general nuclei with 
spin / > 1/2, where Eq. (|10fl does not apply. Then cr(t) 
will be a random walk variable with equal probabilities 
of moving right or left, 

a(t) — 21 A, (21 — 2)A, . . . , — 27A. (26) 

Of course for I > 1/2 we can not represent a(t) as an an- 
alytical function of the number of steps N(t). However, 
for small t (such that Xc(t) ~ 1, or tT~^ <C tA nm <C 1) 
Eq. IL'l'l) will still hold, since a(t) will take at most one 
step (of length 2A). We will use this fact when consid- 
ering GaAs nuclei, which have 1 = 3/2 (see below). 

B. Calculating the nuclear flip-flop rates 

We now show how to calculate the rate T~^\ appear- 
ing in Eq. i|18|). for dipolar coupled spin I nuclei. As 
noted above, we will assume the nuclear dynamics to be 
decoupled from the central spin S through the effective 
Hamiltonian Hj, Eq. JSJ. Furthermore, we will separate 
Ttj in a secular part Ho and flip-flop terms F nm (t) which 
contain an harmonic time dependence: 

H r = Ho + F nm(t), (27) 



to introduce the frequency uj: The rate T~^\ will be ob- 
tained by taking the u> — ► limit. Suppose we have N 
nuclear spins I„ . We will label the eigenstates of Ho by 
the index a — 0, 1, . . . , 2^: Ho\a) = E a \a). Obviously 
we simply have \a) — \m\mi ■ ■ •tun) with rm = ±1/2. 
The transition rate between two of these states induced 
by F nm (t) is then given by Fermi's golden rule to be 

W a T = \bl m MF nm \b)\ 2 [5{E a -E b -u) 

+ 6(E a -E b + , (30) 

with F nm = (I n +/ m _ +J n _i m -|_). The central spin phase 
changes by A nm — \A n ~A m \/2 during one of these 
events. Moreover, we will assume the nuclei in thermal 
equilibrium, each state \a) populated with a Boltzmann 
probability 

P(a) = exp (--^)/£exp (31) 

Then the flip-flop rate for a pair nm becomes 

T m >)=Ep(«) W X (32) 

a , 6 

At zero temperature p(a) will be non-zero only for the 
ground state. Furthermore, if B S> b nm /jil ~ 0.1 G 
this ground state will simply be |0) = \I, I, ■ ■ ■ , I), and 
= f° r au pairs nm, leading to vanishing zero tem- 
perature decoherence. However, nuclear Zeeman ener- 
gies are only ~ 1 mK/Tesla, much lower than typical 
dilution refrigerator temperatures (> 10 mK). Hence a 
reasonable approximation is to assume fc^T 3> hjjB, 
and p(a) — (21 + 1)~ N for all states \a). This in fact 
makes Eq. (|32|l much easier to calculate, particularly 
since we will employ a technique similar to the method 
of moments due to Van VleckiSiS^ Therefore we assume 
p(a) to be the same for all states a, and rewrite Eq. il32*|) 
in the form 

r nmM = 2Kb 2 nmPnm (uj), (33) 

with p nm (u}) playing the role of a density of states given 
by 



Ho = - hi B - ) J «*- 4 £ b nm I nz I mz , (28) 

F nm (t) = b nm (I n+ I m - + I n -I m+ ) cos (cot), (29) 

where we have introduced a fictitious frequency variable 
uj for theoretical convenience - eventually we are inter- 
ested in the uj = limit. We shall see that it is useful 



Pnm(w) = 2(2/ + l)iv El ( a \ F nm\b)\ 2 6 (E a ~ E b ~ Uj) . 

a,b 

(34) 

The invariance of the trace operation allows us to calcu- 
late any moment of p n m(^) exactly. We define the n-th 
moment as 

/>oo / />oo 

(oj n ) = J iv n p nm (uj)duj I J p nrn (uj)duj. (35) 
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As an example of how to calculate these moments, con- 
sider the normalization constant 

/>oo 

A{I) = / Pnm(u)duJ 
JO 

= 4(2J | 1 w E(°i f -i t )^-» 

a,b 

4(2/ + l) N X nmS 

= ^7TT [2/(J + 1) + 1] - (36) 

It is straightforward to prove the following relations: 
(u))=w = 2C(I) Tr {[Ho, F nm ]I n+ I m -} 

= l\A n -A m \, (37) 



C(J) = 



15 (21+1) 



-JV+l 



8 /(/+l)[2/(/+l) + l] : 



w" = - 



C(J) Tr{ [Wo, Fnm] 2 } 



(38) 



(39) 



<(^-^) 2 > = fm+i) E ( & ™- & ™) 2 ' (4°) 



C(I)Tr{[H ,[H ,F nm }] 2 } 



(41) 



((c-ZJ) 4 ) = y/(/+l) 



;[2/(/+i) + i] 



+/(/+!) 



(42) 



Using these expressions we can discuss various possibil- 
ities for the shape of pn m (u)), which is clearly an even 
function of ui with peaks at ±\A n — A m \/2. The most 
common choice in lineshape theory is to assume each 
peak to be a Gaussian or a Lorentzian with cut-off at 
the wings£ To decide between these, we define a dimen- 
sionless parameter 



3<(^-^) 2 ) 2 



i-3A E 

3/ Jr E 



ipni ^mi) 



(43) 



Here we have introduced the occupation probability / 
for a nucleus at site i (sums involving one site and two 



sites are proportional to / and / 2 respectively). For a 
Gaussian function this ratio would be exactly 1. For 
a Lorentzian, it would be equal to ~ a> nm /5 nm 3> 1, 
where a nm is the cut-off and 5 nm the half-width at half 
maximum. Clearly / -C 1 leads to Lorentzian rates. Here 
we will adopt the following approximation: Whenever 
£mn < 10 we will approximate p n m(w) by a Gaussian; if 
£nm > 10 we will use the Lorentzian fit. The Gaussian 
approximation leads to 



Pnm(w) = A(I) 



1 



27TKr, 



exp 



(W - U) 2 



9k 2 



+ exp 



with 



(44) 



(45) 



The desired flip-flop rate is then obtained by setting u> 
in Eqs. lO and (J33JI, 

\A n — A m 



^(U < 10) = 2V2ttA(I)^ exp 



hp 



J nrn 



Fore 



(46) 



> 10 we assume the Lorentzian form 



Pnm,{u) = A(I)- 



7T 

((J — LU + Q',„ 



6>(w 



(w + £J + 



+ 



+ (w - w) 2 
9(w + ZJ 



S'l, 



(47) 



with 9 being the step function. It is important to note 
that this Lorentzian must have a cut-off so that its second 
and higher moments do not diverge. The parameters a nm 
and S nm can be easily related to Eqs. jSP} and (021 : 



2 V / 3V ((^-^) 4 ) 
/2((w-5J) 



4\ ' 



(48) 



(49) 



The Lorentzian flip-flop rate is 

Tn^nm > 10) 



& 2 1 

AA{I)- nm 



0(\A 



1 + [{A n - A m )/(2S nm )Y 
(50) 

Therefore T n i will be suppressed unless \A n — A m \ < 
K-nm , 5 nm . This effect stems from energy conservation: 
For a flip-flop to happen, an amount \A n — A m \/2 of 
energy must be absorbed by the dipolar term in Hq. To 
see this consider two states with 7 = 1/2, 
\a) = | + l/2,-l/2,m 3 ,...,TOjv}, 

1 6) = | — 1/2, +1/2, 777-3, ■ ■ ■ , Trijf). Their energy difference 
is 



E a - E b = \ {A x -A 2 )-4J2 (hi ~ h) : 



(51) 
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If we calculate the average square of this expression using 
(to?) = 1/4, (to») = wc get 



S b ) 2 ) = ( -^-^ ) • 1 X (6« - 6 2 «) 2 , (52) 



showing that an adjustment of the spin system can sup- 
ply energies up to ~ K nm to compensate for the cen- 
tral spin field. We believe that the above expressions for 
are considerably more precise than the phenomeno- 
logical rates used in the earlier literature without any 
derivations £iL2i 



III. APPLICATIONS 

A. Electron spin of a phosphorus donor in silicon 

We now turn to applications of our theory for systems 
of interest to quantum computation in spin qubit-based 
semiconductor architectures. We start by considering 
the electron spin of a shallow donor in silicon. In nat- 
ural samples, 95.33% of silicon atoms have no nuclear 
magnetic moment^ Those are the 28 Si isotopes. Spec- 
tral diffusion is then caused by the remaining fraction 
fnat = 0.0467 of 29 Si isotopes which are spin 1/2 nu- 
clei with gyromagnetic ratio jf 1 = 5.31 x 10 3 (s G) _1 . 
These nuclei will produce a hyperfine field on the electron 
donor impurity given by Eq. JSJ , which is proportional to 
the electron's probability density at the nuclear site R„, 
|\&(Rn)| 2 . For this state we assume the Kohn-Luttinger 
wave functioni 3 - 2 ^ 



*( r ) = -iE^'( r K-( r ) e<krr ' 

V .7 = 1 



(53) 



2lT - 

0.85 kj, kj € {x, — x, y, — y, z, — z} ,(54) 

asi 



exp 



+ 



V 2 +z 2 



(nb) 2 (na) 



y/Tr(na) 2 (nb) 



(55) 



with the appropriate corresponding envelope functions 
Fj [Eq. 1)55(1 ] with anisotropics in the y and z directions. 
Here n = (0.029eV r /£' i ) 1 / 2 with Ei being the ionization 
energy of the impurity (Ei = 0.044 eV for the phosphorus 
impurity, hence n = 0.81 in our case), agi = 5.43 A the 
lattice parameter for Si, a — 25.09 A and b — 14.43 
A characteristic lengths for Si hydrogenic impurities. 13 
Moreover, we will use experimentally measured values 
for the charge density on each Si lattice sitejSSiSI 



Hence the hyperfine interaction is given by 

A n = ^TsSf hv[Fi{^n)cos(k X n ) 
y 

+F 3 (R„) cos (k Y n ) + F 5 (R„) cos (k Z n )} 
1 - 3 cos 2 9 



-7s 7/ n>- 



|R„| 3 



|R„| - na), 



(57) 



with k = (0.85)27r/asi, and 7 | l = 1.76 x 10 7 (s G)" 1 the 
gyromagnetic ratio for the electron donor. It is instruc- 
tive to check the experimental validity of Eq. (|57|l by 
calculating the inhomogeneous line-width (~ l/7| l T 2 *). 
A simple statistical theory applied to Eq. (JJJ leads toSi 



{^h s s l -By 



1 . 



(58) 



For / = 0.0467 our calculated root mean square line- 
width is equal to 0.89 G. On the other hand, an ESR 
scan leads to 2.5 G/2V2 In 2 = 1.06 G; 1 - 3 Therefore our 
model is able to explain 84% of the experimental hyper- 
fine line- width [the residual dipolar term in Eq. I|57|) only 
contributes ~ 0.1% to this line- width]. 

Before discussing our spin echo decay results we should 
mention how the 29 Si fraction / enters our calculations. 
For example, Eq. I|14f> becomes 



v(2r) = n [ (2t)] 



(59) 



since the probability of a pair n, to to be 29 Si is / 2 . Also 
single and double sums in Eqs. I|40|) and (|42|l are propor- 
tional to / and / 2 (these sums are calculated assuming 
the pair n, m is already occupied; hence the question is 
whether or not site i ^ n, m contains a 29 Si). To perform 
our numerical calculations we take the natural logarithm 
of Eq. (|59|l . Then we take advantage of the fact that 
b nm oc [Eq. ©] decays fast as a function of the 

inter- nuclear distance R nm - Therefore we achieve faster 
convergence by summing over lattice sites n together with 
some of its close neighbors to. Both the range of sites n 
and the number of neighbors included in the sum were 
increased systematically to ensure proper convergence. 
From the convergence check we concluded that summing 
within eight characteristic lengths of the wave function 
(~ 40asi for the Si:P case) plus including up to six near- 
est neighbor shells gave excellent convergence. 

Fig. 1 shows the echo signal [Eq. (|18fl ] as a func- 
tion of time for natural Si and B \\ [111]. Clearly the 
echoes occur at t = 2r. The echo envelope, defined by 
the maximum of each peak, should be compared with the 
empirical experimental fit of Chiba and Hirai, 17 



VExpi^r) = exp 



— V(n£-Y 
0.6ms J \ 0.4ms / 



(60) 



\uj(R r , 



Tj w 186. 



(56) 



Clearly our theory is able to explain quite successfully 
the exp (— r 3 ) decay, but our exponential tail is twenty 
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Time [ms] 

FIG. 1: Calculated echo peaks VE(t) and echo envelope 
ve(2t) as a function of time for a single donor electron spin 
in Si:P. We assumed the Si lattice had natural abundance of 
29 Si (/ = 0.0467) and B || [111]. We also show a fit for our 
echo envelope that compares well with the experimental fit of 
Ref.[l3 [see Eq. iSSH . 



times smaller than the measured value. Therefore our 
theory suggests this extra exponential decay is coming 
from other decohcrcnce mechanisms, perhaps related to 
imperfections of the ESR pulses [Recent experimental 
data— suggest dipolar scattering between donor electrons 
is responsible for the extra exponential decay seen in Eq. 
(|60|l ]. Nevertheless, our Tm — 0.64 ms is ~ 2 times larger 
than the measured value of 0.3 ms, a quite good agree- 
ment in view of former spin relaxation calculations— In- 
cluding the residual dipolar term in A n [see Eq. (|57|l ] 
changes Tm by less than 1%, showing that here SD is 
dominated by flip-flopping nuclei inside the electron's 
wave function [if we set the hyperfine term to zero in 
Eq. (52J, we get T M = 1.7 - 1.9 ms for < r < no]. 
It is interesting to note that former SD estimates^iii ne- 
glected the hyperfine contribution, arguing that nuclei 
inside the electron's wave function could not flip-flop due 
to the large A nm they would induce. Surprisingly, our 
calculation shows that this coupling is very important, 
although estimates using only the dipolar term would 
still lead to reasonable results. The hyperfine term be- 
comes more and more important as the wave function 
size increases. Fig. 2 shows the behavior of — 1iiu.e(2t) 
for B || [111] and a few values of /. As 2r is increased, 
the echo envelope changes from a "Gaussian spectral dif- 
fusion" regimeii [~ r 3 , well described by Eq. (|22|) ] to 
-lav — r [see Eq. (JT^ and for A nm r > 1]. The 
"experimental window" 0.1 < — lav < 3 falls within 
this crossover for many values of /, including natural 
isotopic abundance. The exponential behavior obtained 
for long 2r can be understood by looking at the corre- 
lation function [Eq. 124fl ] which also appears in Fig. 2 
(right scale) . Notice that this correlation function decays 
smoothly over 3 time decades, but its behavior is clearly 
non-exponential (see Fig. 9), making evident the dif- 




10" 5 10" 4 10" 3 10" 2 10" 1 10° 
2x [seconds] 

FIG. 2: The left scale shows the logarithm of the -k/2 — tt echo 
envelope as a function of 2r for electron spins of Si:P. For all 
29 Si isotopic fractions / this function undergoes a crossover 
from Gaussian spectral diffusion behavior [v ~ exp(— fcr 3 )] 
to motional narrowing [v ~ exp (— fcr)]. The correlation func- 
tion, Xc(2t) is shown on the right scale; motional narrow- 
ing behavior corresponds to small values of \ a . This corre- 
lation function shows non-exponential decay, quite distinct 
from former theories^*— *— (see also Fig. 9 below). Here we 
have B || [111]. 

ference between our current complete theory and former 
heuristic ones based on the two parameter Gaussian con- 
ditional probabilitypiiiSiii The small correlation between 
Zeeman frequencies before and after the application of 
the 7r pulse explains this motional narrowing behavior 
for v(2t) (see discussion in the end of section II- A). Note 
that Fig. 2 suggests it might be more appropriate to 
fit — \av(2r) to a function which approaches r at long 
times and t 3 at short time scales. Such a fit could pos- 
sibly yield a better description for experimental data as 
long as other contributing mechanisms (electron-electron 
dipolar scattering^ for example) are properly subtracted. 

It is interesting to see if the Si six-fold degeneracy em- 
bedded into the Kohn-Luttinger state [Eq. Qfnljl] has a 
strong effect on Tm- We obtain Tm ~ 0.3 ms at natu- 
ral abundance using a hydrogenic state with Bohr radius 
~ 20 A. Therefore we can conclude that the strong oscil- 
lations of A n reduce T~^ by a small amount, hence in- 
creasing Tm (as A nm is larger the energy cost for nuclear 
flip-flops is higher). This gives an idea about the effect 
of the electron's field on the 29 Si nuclei. Fig. 3 shows 
Tm as a function of /. Only one purified sample with 
/ = (0.12 ± 0.08)% was studied experimentally^ lead- 
ing to Tm — 0.52 ms. Note that this value is significantly 
different than what we have in Fig. 3 (~ 0.1 s) simply be- 
cause at such low values of / the concentration of donors 
(4 x 10 16 cm~ 3 ) lead to strong electron-electron dipolar 
scattering (see Fig. 2 of Ref. @). It would be quite inter- 
esting to see spin echo measurements for different values 
of /, but as the isotopic purity increases one will have to 
decrease donor concentration significantly to ensure that 
dipolar scattering is not prevailing (for that purpose Tm 
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FIG. 3: Decoherence time (Tm) of the Si:P electron spin as 
a function of 29 Si content /. The straight line resembles a 
simple theory derived previously- where / only altered the 
flip-flop pair probability. Our new theory deviates from this 
behavior since / also enters the moments [Eqs. (M0> and l|42|l ]. 
The change in slope seen in Tm is explained by the transition 
from Gaussian to Lorentzian flip-flop expressions for T~^ [See 
Eqs. (gSJ and CT l. 



should be independent of donor concentration, showing 
saturation for low P concentration^). 

An interesting confirmation of our theory would be to 
measure Tm as a function of the B field tilting angle 9. 
Fig. 4 shows this angular dependence in samples with 
natural isotopic abundance. 8 = 0° means B || [001], 
while 9 = 90° B || [110] (unfortunately Ref. IftI measured 
Tm only for B || [111]). Preliminary experimental data 
shows excellent qualitative agreement with Fig. Inci- 
dentally, our theory does not depend on the B field inten- 
sity, except for B > 10 T: Then k B T > fry/ 5 does not 
hold, and also the magnetic length Ib = (He/ eB) 1 / 2 < 20 
A, and B is already deforming the electron's wave func- 
tion. The latter effect does in fact appear in the quantum 
dot case, since its radius depends crucially on Ib [see Eq. 



B. Nuclear spin of a phosphorus donor in silicon 

31 P nuclear spins [7^ = 1.08 x 10 4 (s G) _1 ] are promis- 
ing qubit candidates when implanted in a Si matrix^ Its 
interaction with the lattice is rather weak, leading to a 
measured T\ ~ 5 hours at B = 0.8 Tesla and T = 1.25 
K>i4 Here we calculate the contribution of 29 Si spectral 
diffusion on the decoherence time Tm of an isolated 31 P 
nucleus. There is an important difference between this 
calculation and the one related to the Si:P electron spin 
above. Even though the flip-flop rate of a pair n, m of 29 Si 
(Tnm) i s determined in exactly the same way as above 
[Eqs. (O and J3UJ) with A n given by Eq. (JSJ)], the 



FIG. 4: Behavior of Si:P electron spin Tm as a function 
of B field tilting angle with respect to the crystal lattice for 
samples with natural abundance of 29 Si. To our knowledge, 
this dependence has never been probed experimentally, and 
would be an interesting test for the accuracy of our theory. 



corresponding phase change on the 31 P spin is given by 



\A' n — -Am 1/2) 



1-3 COS 2 ( 



(61) 
(62) 



where A' n is the dipolar interaction between a 29 Si nu- 
cleus located at R n and the 31 P located at the origin (in 
principle we should include the effect of the 31 P dipo- 
lar field on the flip-flop rates T~ 7 \, but we determined 
it to be negligible on the Tm calculation both in this 
section and the preceding one). It is easy to see that 
this phase change A nm is in most cases ~ 10 3 smaller 
than the corresponding one for the electron. Therefore 
the motional narrowing condition 3> A nm will be 
satisfied for many more pairs than in the last section. 
Hence we should expect v ~ exp (— fcr) for a wide range 
of parameters [see Eq. (1191) ] . 

Fig. 5 shows the calculated shape of the echo envelope 
as a function of It. The qualitative behavior is similar 
to Fig. 2 above, except for the rather abrupt crossover 
from Gaussian SD to motional narrowing behavior. For 
/ < 0.0467 the observed decay is indeed v ~ exp (— kr) 
as predicted above. This behavior can be understood by 
looking at the correlation function \c (right scale on Fig. 
5) which goes to zero in only one time decade, evidencing 
the abrupt appearance of the motional narrowing regime. 
Unfortunately there are no experimental data available 
to verify this qualitative result. Fig. 6 depicts the de- 
pendence of Tm as a function of /. 

The dependence on the tilting angle 9 (Fig. 7) is sig- 
nificantly different than the one in Fig. 4. This can 
be explained by the fact that while in the former case 
first nearest-neighbor flip-flop dominated (about 92% of 
T^ 1 ), here first nearest neighbor amounts only to 3% of 
the rate, which is dominated by second nearest neighbors 
(~ 90% of Tm 1 )- This is why we see more oscillations in 
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FIG. 5: Logarithm of 7r/2 — ir echo envelope as a function 
of echo time 2r for 31 P nuclear spins in Si:P. This function 
undergoes an abrupt crossover from Gaussian spectral diffu- 
sion behavior to motional narrowing, quite different from the 
smooth transition seen on Fig. 2. On the right scale we plot- 
ted the correlation function Xc, which shows a sharper tran- 
sition to motional narrowing as compared to Fig. 2 (here \ c 
changes over one time decade, as opposed to three in Fig. 2). 
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FIG. 6: Depicts Tm versus Si content / for P nuclear 
spins. Similar to Fig. 3, the change in slope around / ~ 0.01 
occurs due to the crossover to Lorentzian T^ m [Eq. 



Fig. 6 than in Fig. 3: Second nearest neighbors are 12 in 
number, and have a more intricate lattice configuration 
than the 4 first nearest neighbors. Quite interestingly, in 
Fig. 7, Tm is maximized when B points in the [110] di- 
rection, and displays a sharp peak in the [111] direction, 
while Fig. 4 shows Tm minimized when B points in both 
directions. 



C. Electron spin in a gallium arsenide quantum dot 

Recent interest in the spin properties of a single elec- 
tron GaAs QD is motivated by its potential use as a 
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FIG. 7: Spectral diffusion decay time Tm versus B field tilt- 
ing angle for 31 P nuclei. This calculation assumes natural 
abundance of 29 Si, / = 0.0467. 



qubitiSS One advantage over Si:P is the non-degenerate 
conduction band, with the minimum at the T point. 
Therefore the exchange interaction of a double dot is a 
smooth function of the inter-dot barrier and distance, 
quite different from two Si:P donors, where this ex- 
change may oscillate dramatically as a function of donor 
separation^ 

The GaAs lattice has a Zincblende structure with 50% 
of 75 As, 30.2% of 69 Ga, and 19.8% of 71 Ga, 2& These nu- 
clei have spin / = 3/2, but here we present calculations 
using the spin 1/2 stochastic theory together with spin 
3/2 flip-flop rates. Since the correlation function is very 
close to 1 in the neighborhood of our calculated Tm, the 
electron is well into the slow SD regime, making this ap- 
proximation very accurate for our purpose. The electron 
hyperfine interaction with these nuclei will lead to an 
inhomogeneous linewidth of about 50 G for small dots 
(with Fock-Darwin radius £ ~ 20 nm; donor impurities 
in bulk GaAs have even higher broadening 31 ). If one has 
an ensemble of these dots, the decoherence time Tm can 
be measured by using a 7r/2 — ir pulse. However, en- 
sembles of dots always contain size distribution. Since 
Tm will be quite sensitive to the radius I (see below), 
it might be more appropriate to measure decoherence by 
applying an ESR field to a single dot and then measuring 
the signal using transport experiments, even though such 
experiments would only lead to a lower bound on Tm? 2 

The energy eigenstates of a quantum dot in the pres- 
ence of spin-orbit coupling are a mixture of spin states 
up and down. Therefore if the electron spin is up, it will 
flip within a time T± with the corresponding emission of 
a phonon. However, for the spin to flip a virtual tran- 
sition to an excited orbital state has to happen, since 
spins do not couple directly to the phonon strain field. 
The result is a strong sensitivity on dot size and applied 
field B, Tf 1 cx £ 8 B 5 . For B ~ 1 Tesla and I ~ 30 nm 
a recent theory leads to T\ ~ 1 ms, showing that for the 
small dots in a quantum computer architecture spin-flip 
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scattering is strongly suppressed, 20 Recently we showed 
that the dominant decoherence mechanism in these small 
dots (£ < 50 nm) is nuclear spectral diffusion. 6 Here, our 
detailed calculation of this effect confirms the accuracy 
of the simple theory presented earlier^ justified by Eq. 
(|22|l and the fact that the correlation function [Eq. I)24|l] 
is still quite close to 1 in the neighborhood of Tm- For 
\z\ < z q/2, the quantum dot wave function can be simply 
approximated as 



W r = \ — cos — z\ _ ; — - 



exp 



x 2 + y 2 \ 
2£ 2 {B) ) 



(63) 

and we assume ^ = for z > Zq/2. This state is a reason- 
able description for the lowest orbital of a quantum well 
of thickness zq, with electrostatic lateral parabolic con- 
finement with radius do- The Fock-Darwin radius t(B) 
includes the additional B field confinement, 



1(B) = 



Wo 




(64) 
(65) 



Then the hyperfine coupling A n for a nucleus located 
at a coordinate (X n ,Y n , Z n ) from the center of the dot 
becomes 



An = l 4 1Sl J h f? aAS d{I)cOS 2 [-Z n 



3 £ 2 (B)z Q 
x exp 



X: '"" h: *0(za/2-\Z n \) (66) 



£ 2 {B) 
1-3 cos 2 e n 

irj 5 



9(X 2 +Y 2 -£ 2 (B)).(67) 



Here a GaAs = 5.65 A, 7s = -3.86 x 10 6 (s G)' 1 (assum- 
ing g — —0.44 independent of I for the dot electronic), 
7/ = 4.58, 8.16, 6.42 x 10 3 (s G)" 1 and charge densities 
d{I) = 9.8 , 5.8 , 5.8 x 10 25 cm- 3 for 75 As, 71 Ga, and 69 Ga 
respectively^ 3 " The residual dipolar coupling again leads 
to a very small correction (< 1% of T M ). 

Fig. 8 shows the behavior of Tm as a function of Fock- 
Darwin radius i for three different quantum well thick- 
nesses z . For each z , Tm displays a minimum as a func- 
tion of £. This arises from two competing effects: If we 
decrease the wave function size, fewer nuclei contribute 
to SD, hence Tm — f oo as I — > 0. However when we 
increase £ the wave function flattens, making two close 
pairs n, m produce similar hyperfine fields. Since it is 
the difference \A n — A m \ which causes phase fluctuation, 
Tm ^ oo as £ — > oo. Therefore Tm must have a mini- 
mum as a function of the Fock-Darwin radius £. At this 
minimum, Tm ~ 10 /is for all zq. For all QD sizes consid- 
ered here the SD decay is found to be ~ exp (— t 3 ), the 
correlation function being close to 1 in the neighborhood 
of Tm (Fig. 9). The dependence with 9 is similar to Fig. 
4. 



40 



^30 



20 



1 / 

/ 

/ 


/ Quantum Well 




thickness, z 








— ■ — 1 nm 




—A— 20 nrn^-- 1 







50 100 150 

Fock-Darwin radius [nm] 



200 



FIG. 8: Tm versus Fock-Darwin radius I for various quantum 
well thicknesses, zo = 5, 10, 20 nm. Decoherence achieves a 
maximum as a function of I. Here B \\ [111]. 



IV. CONCLUSION 

We have developed a detailed quantitative theory for 
nuclear spectral diffusion of localized spins in semicon- 
ductors. By treating each nuclear pair independently, we 
are able to show that the echo signal arises from the in- 
terplay between fast and slow nuclear flip-flops. The echo 
envelope undergoes a smooth transition from a Gaussian 
SD decay to an exponential motional narrowing when the 
nuclear bath loses correlation over time, this transition 
being well described by an appropriate correlation func- 
tion. The Lorcntzian approximation gives a good descrip- 
tion of the intermediate crossover regime, and our theory 
gives a microscopic justification for the use of these phe- 
nomenological conditional probabilities. This behavior is 
quite general and should be observed in any decoherence 
mechanism where qubit phase fluctuation takes place. 
We apply our theory to three physical systems proposed 
as QC architectures, showing that SD should not be a 
decisive constraint in their development (Tm > 10 /is; 
therefore Tm/tj 3> 10 6 for all architectures, where tj is 
the time scale for the exchange gate 6 ). Our calculated 
Si:P electron donor spin results agree well with existing 
experimental data for natural abundance of 29 Si, while 
Tm increases very fast as these isotopes are removed from 
the lattice. In addition our calculation shows that the 
most important contribution to SD comes from nuclei 
inside the electron's wave function spread, as opposed to 
electrons far away as was suggested before H The 31 P 
nuclear spin is found to be in the motional narrowing 
regime, weakly affected by 29 Si. Our GaAs quantum dot 
calculations confirm our earlier estimates based on a sim- 
pler theoryi 6 - Although there are no Ga or As / = iso- 
topes, one way to reduce SD is to suppress flip-flop events 
by nuclear polarization. The main difference between T\ 
and T2 samples is the presence of several relaxation rates. 
This feature is evident in Fig. 9, where we show corre- 
lation functions for many cases treated here. Increasing 
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away from 2r as long as a single pair undergoes slow SD. 
A global treatment of all pairs such as the one achieved in 
the two parameter Gaussian theory leads to an echo peak 
formed between r and 2tiAI Nevertheless to our knowl- 
edge this effect does not seem to be present in 7r/2 — tt 
echo experimentsjSiii our approximation being appropri- 
ate in this respect. Another approximation that seems 
reasonable is to neglect back action. We assume that the 
nuclei evolve under a static central spin field [Eq. (JSJ], 
unchanged by its effect on this spin. 
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FIG. 9: Non-exponential decay of the correlation function \c 
as a function of time [see Eq. 123H 1. Initially \c shows strong 
exponential decay, dominated by the fastest rate. Then slower 
processes take over. For Si:P we assume natural 29 Si abun- 
dance, while B || [111]. GaAs-QD's show a clear tendency to 
motional narrowing as £ (Fock-Darwin radius) increases. 



QD radius I increases Tm and reduces Xc{Tm), pushing 
large QD's to the motional narrowing regime. This result 
is clear evidence that SD does not affect delocalized states 
(such as conduction electrons), since all nuclear pairs will 
be in the motional narrowing regime [Eq. (|19fl ] . 

We now briefly discuss the approximations assumed in 
our work. Certainly the most important one is that we 
treat nuclear pairs as independent Markovian processes, 
making random flip-flops with a Poisson distribution. We 
believe this approximation is reasonable for a nuclear 
bath at temperatures much higher than its dipolar order- 
ing critical temperature (~ hb nm /kB ~ 10 -9 K), since in 
this case the spin correlation length can be assumed of 
the order of one lattice parameter. Of course this con- 
dition will always hold for dilution refrigerators (which 
operate at millikelvin temperatures), and specially for 
the experiment discussed in this paper (T — 1.6 K)ii At 
high spin temperature the nuclear system undergoes fre- 
quent transitions through its (21 + 1) N available states. 
An important assumption implicit in our model is that 
only a subset of these transitions, namely flip-flops be- 
tween close nuclear pairs [described by a nm , Eq. J5J] are 
responsible for the SD effect at external magnetic fields 
B S> 100 G. Four spin processes stem from higher order 
perturbations in Eq. Q2'jp. and can be neglected when 
b nm < A„ m , which is the case for the relevant nuclear 
pairs [otherwise, the condition b nm ^> A„ m immediately 
implies motional narrowing [see Eq. I|19(l ] and this pair 
gives a negligible contribution to the SD rate]. Prelim- 
inary experimental results^ on the orientation depen- 
dence of Tm show excellent qualitative agreement with 
Fig. 4, suggesting the pair flip-flop picture is quite appro- 
priate. One limitation of the uncorrelated pair approx- 
imation is that the echo peaks always occur at t = 2r. 
This happens because the product [Eq. I|14|) ] will be zero 



Finally, we comment on the relationship between 
our current work and a number of other recent pub- 
lications in the literature dealing with spin relaxation 
in the context of semiconductor quantum computer 
architecturesi^SiSii^! We emphasize that our theory deals 
exclusively with electron spin decoherence in the absence 
of any phonon effects. Our work provides a theory for 
nuclear-induced spectral diffusion of (localized) electron 
spin in semiconductors, and as such we deal with elec- 
tronic spin decoherence arising from the fluctuations in- 
duced by electron-nucleus hyperfine interactions caused 
by nuclear spin dipolar flip-flops. Since the flip-flops con- 
serve nuclear energy, no phonons are required for the de- 
coherence studied in this work. 



Our work is the first comprehensive theory for nuclear 
induced spectral diffusion in the spin decoherence of lo- 
calized electrons in semiconductors - we do not invoke the 
empirical approximations and the heuristic arguments 
which were used to obtain earlier expressions for spec- 
tral diffusion existing in the literaturepii*ii24 As such 
our results are applicable and relevant not only to consid- 
erations involving solid state spin-qubit-based quantum 
computer architectures but also to all problems involv- 
ing spin decoherence due to electron-nucleus hyperfine 
coupling where phonon effects are negligible (i.e. at low 
temperatures). In particular, our results apply to spin 
echo measurements in semiconductors at low tempera- 
tures, and it is therefore gratifying that we have been able 
to quantitatively explain hitherto unexplained Si spectral 
diffusion results of Ref. 17 dating back thirty years. We 
emphasize that our theory is still approximate since we 
are forced to make a number of approximations, the most 
important one being the assumption of uncorrelated flip- 
flops among spin pairs. Although we believe this assump- 
tion of uncorrelated flip-flops to be well- valid at "high" 
nuclear spin temperature (3> nK), it is still worthwhile 
to consider further improvement of our theory by taking 
into account the full non-Markovian nature of the spin 
flip-flop processes in future work. The authors acknowl- 
edge discussions with S.E. Barrett, A. Kaminski, S.A. 
Lyon, J. Fabian, and P. Zoller. This work is supported 
by ARDA, LPS, US-ONR, and NSF. 
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APPENDIX A: STOCHASTIC THEORY FOR A 
SINGLE FLIP-FLOPPING PAIR 

The problem is to evaluate the average [Eq. (|15Jl ] for 
a single nuclear pair n, m 



v(t) 



cos 



A / S(t')(-l) N ^dt' 
o 



(Al) 



with N(t) a Poisson random variable with parameter t/T 
(for simplicity we dropped the subscript nm from v nm , 
A nm , T nm ). Expanding the cosine and rearranging the 
product of integrals we get 

v(t) 



■ — ' nZ 

V(-l) fe A 2fc / dt 2k S(t 2k ) 

rt2k 

X / dt 2 k-lS{t 2 k-l) ■ ■ ■ 
Jo 

Jo 



(A2) 



with £ = N(ti) + ■ ■ ■ + N(t 2k ). Using the inequality 
t > t 2k > t 2 k-\ > • • • > t\ > together with the fact 
that a sum of two Poisson variables with parameters 4 j /T 
and tj/T equals another Poisson with (4j + tj)/T, we get 

«e = N(h k -t2k-l) + N(t2k-2-t 2 k-3) + --- 

+N(t 2 -h) + 2[N{t2k-i) 

+N(t 2k -3) + --- + N(h)], (A3) 

((-l) 5 ) = expj - ~ - tak-i) + (*2fc-a - *2*-s) 



T 



+ -- + (ta-*i)] 



(A4) 



since the number of flip-flops N(t 2k — <2fc-i) are assumed 
independent random variables for non-overlapping time 
intervals £2/0 — t 2k -i (Markovian approximation). Using 
Eqs. I|A2() and (|A4|) we rewrite the echo decay in the 
form 



;(*) = £(-l)*A 2 *t»*(t) 



(A5) 



fe=0 



with v 2k (t) satisfying the integral recurrence relation 

ri 



V 2 k(t) 



t 



X 



di 2 fe-iexp ( 2 t2 ^ 1 )S , (t 2 /c-i) 



T 



X«2fc-2(*2fe-l), 



(A6) 



and i>o(t) = 1. We can transform Eq. (|A6|) into an 
algebraic recurrence relation by using the properties of 
Laplace transforms^ 



v 2k (p) 



-pt 



v 2k (t)dt. 



(A7) 



The free induction signal, vp(t) [Eq. I|17f) ] is obtained by 
setting 5(4) — 1 for all times. Then Eq. (jA6(l becomes 



^(P) = p( p + 2T -i) ^-^)> 



(A8) 



and the Laplace transform of VF(t) is easily calculated to 
be 



vf(p) 



P+2T- 



p( P + 2T- 1 ) + A 2 



(A9) 



The inverse transform of this expression can be obtained 
by expanding in partial fractions^ 



vp(t) = exp 



T 



RT 



sinh (Rt) + cosh (Rt) 



(A10) 

where R 2 = T 2 — A 2 . We now turn to the tt/2 — tt echo 
w.e(4), which is obtained by setting the echo function to 

S(t) = 1 - 20(4 - t). (All) 

Then Eq. (|A6|) becomes 

^2fc0) = p (p + 2 T -i) «2fe-2(p) - 2exp(-pr) 

xexp(-2TT" 1 ) / eft' exp (24'T -1 ) 
Jo 

xv 2k - 2 (t') (A12) 

and after summing the series we get 
P + 2T- 1 



v(p) 



A 2 +p{p + 2T- 1 ) 
2A 2 exp (-pr) 



A 2 + p(p + 2T" 



k/(t), (A13) 



/(r) = exp(-2rT- 1 ) 

x / d4'cxp(24'T" 1 )t; J .(4')- (A14) 
Jo 

Using Eq. l|AT0)l . the inverse transform becomes 

2A 2 

u(4) = v F (4)+0(4-r)— 2- exp(-4T- 1 ) 

x sinh (Rt) sinh [i?(4 - t)] (A15) 
= 9(T-t)v F (t) + 6(t-T)v E (t), (A16) 



which after rearrangement leads to Eq. Ijl8|l for t > r, 

u B (4) = i?~ 2 exp (-tT^ 1 )!^ 2 cosh (Rt) 
+i?T~ 1 sinh (i?4) 

-A 2 cosh[-R(4-2r)]l. (A17) 
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APPENDIX B: EQUIVALENCE OF ENSEMBLE 
tt/2 - 7T ECHO ENVELOPE AND AVERAGE 
SINGLE SPIN DYNAMICS 

Eq. (|15fl was derived by averaging a nm (t = 0), which 
clearly applies only to an ensemble of spins. Therefore 
an interesting question is how the single spin off diago- 
nal density matrix element (S±) [analogous to Eq. (|12fl . 
the in-plane magnetization in the S z basis but without 
Cmn(0) average] behaves under nuclear dynamics aver- 
age only. The result is that the modulus squared of this 
quantity is exactly equal to the tt/2 — tt echo envelope, 
as we show below. The evolution under Eq. (JSJ is given 
by 



where Vnm(t) is the free induction decay derived above 
[Eq. I|17|l]. Clearly the difference between a single spin 
and an ensemble is the presence of the complex term in 
Eq. i|B2l) , which leads to strong interference effects when 
the product over pairs is taken. To see this, we calculate 



,(*■) 



(*)] 8 



R 2 



■ exp 



2t 



sinh 2 (R nm t) 



,( E ) 



(2t). 



(B3) 



(Sx(t)) = J] ( ex P lAnm J Q (-l) JV " m( *' ) rf* / 



, (Bl) 



where N nm (t) are independent Poisson random variables 
representing the nuclear dynamics, on which the average 
is taken. Using the same methods of Appendix A, it is 
easy to show that the argument of the product (|B1(I is 
given by 

SLM = (cos A nm j\-l) N ^dt' \ 
A 



sin 



' ' / ( 1 ) 

o 



vi F 2(t) 



+i— — exp 



Tn 



sinh (Rnmt), (B2) 



Therefore the effect of this complex part is to enhance 
the coherence of the single spin, making it exactly equal 
to the tt/2 — 7T echo envelope, 



\(S ± (t))\ 2 =ve(2t). 



(B4) 



Notice that if we performed the product over pairs with- 
out this complex part, S± would decay similarly to free 
induction, vp(t). Many authors define a coherence time 
t c equal to the 1 /e decay of the modulus of the off diago- 
nal density matrix. If the echo is dominated by Gaussian 
SD, we simply have r c = T^f/v^4. 
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